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We present the latest results on event-by-event fluctuations of charged particle multi- 
plicity in Au-|-Au and Cu+Cu collisions at ^/sNN = 200 GeV and 62.4 GeV measured by 
^ ■ PHENIX experiment at RHIC. The two particle correlation length to discuss order param- 
^ . eters, which can be extracted from scale dependence of the fluctuations, is supposed to be 
'nI" ! sensitive to the critical points of QCD phase transition. The obtained correlation length 
^ I with respect to the pseudo rapidity gap indicates the power law behavior as a function of 
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the number of participant nucleons for Au+Au collisions at ^/snn = 200 GeV. 
1. Introduction 



Event-by-event fluctuations on charged particle multiplicity and transverse momentum 
O I have been studied to understand the properties of QCD phase in heavy-ion collisions. 
^ I Several candidates of order parameters in QCD phase, for example correlation length [ 
H] and specific heat [12 E], can be extracted from these observables. It enables us to 
characterize the critical point of the phase transition by surveying the behavior of order 
parameters as a function of temperature or energy density of the system. This is one of 
I the basic methods in the condensed matter physics. We have measured charged particle 
multiplicity fluctuations in Au+Au and Cu+Cu collisions at ^/snn = 200 GeV and 62.4 
GeV by using the PHENIX central arm spectrometer located at mid rapidity region 
(|?7| < 0.35) [13]. The extraction of two particle correlation length has been performed 
from the multiplicity fluctuations measured in various sizes of pseudo rapidity gaps by 
using an integral of two particle correlation function. 

2. Charged particle multiplicity distribution 

In an approximate sense, multiplicity distributions of charged particles produced in sev- 
eral high energy collisions such as e~^e~ , p+p collisions [ 5J and nucleus- nucleus collisions 
[in] are well agree with the negative binomial distributions (NBD) defined as 

p(k) _ r(n + fc) f fx/k \ 1 

where /i indicates the average multiplicity. The finite k value in Eq. describes the 
deviation from the Poisson distribution, because NBD corresponds to the Poisson distri- 
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bution in the limit of infinity for the parameter k. In this respect, NBD k parameter 
represents the muhiphcity fluctuation. 

We have confirmed, that the multiplicity distributions measured by PHENIX agree 
with the NBD for all collision centralities. Fig. ^ shows the multiplicity distributions 
measured in the subdivided acceptances of the pseudo rapidity and the NBD fits to the 
data in — 10% centrality. The collision centralities are defined by the forward detectors 
in PHENIX [ 4,. These are independent subsystems from the central arm detectors. 
The geometrical acceptance is divided into 1/8 through 8/8 of the total acceptance of 
At] = 0.7. For the azimuthal angle, the acceptance is selected as A0 < l/27r. Events 
collected without magnetic field are used to enhance charged particles with low momenta. 



Figure 1. Multiplicity distributions for 
the various pseudo rapidity gaps mea- 
sured by — 10% centrality in Au+Au 
at y/SNN = 200 GeV. The horizontal axis 
is normalized by the average multiplicity. 
Pseudo rapidity gaps are divided into sev- 
eral bins as indicated. Each distribution is 
scaled by the factor indicated. Solid lines 
are fit results with NBD. 

n/<n> 

Fig. 121 shows the measured NBD k parameters as a function of average multiplicity 
and the number of participant nucleons in Au+Au and Cu+Cu collisions at ^snn = 200 
GeV and 62.4 GeV. These are measured with the magnetic field. One can find that the 
NBD k parameters are not scaled by the average multiplicity but scaled by the number 
of participant nucleons for the different energies in Au+Au collisions. 
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Figure 2. Measured NBD k parameters as a function of the average multiplicity (left) 
and number of participant nucleons (right) at the transverse momentum range from 0.2 
to 2.0 GeV/c. Markers correspond to the collision systems and energies as indicated. 
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Figure 3. Corrected NBD k parameters as a function of pseudo rapidity gap measured in 
0% through 70% (left) and 5% through 65% (right) colhsion centrahty with 10% interval 
respectively. Error bars indicate statistical errors and boxes indicate systematic errors. 
Solid lines are fits to Eq. (0). 



3. Extraction of correlation length 

The correlation length ^ in the unit of pseudo rapidity can be extracted from the St] 
dependence of NBD k parameters by assuming a correlation function, 

^2(1/1, 1/2) = p2{yuy2) -pi(i/i)pi(i/2), (2) 

where p2 and pi are the inclusive two particle and the single particle density respectively. 
We introduced a normalized correlation function defined and parametrized as 

R2{yu y2) ^ = ae-\y^-y^\'^ + b, (3) 

pi{yi)pi{y2) 

where ^ is the two particle correlation length and a is the strength of two particle cor- 
relations. Eq. (jS)) dose not have a power term with respect to the pseudo rapidity gap 
size, which is generally introduced to discuss the correlation length around the phase 
transitions, but it is a reasonable formulation for the one dimensional observable [ [7] . 
The two particle correlation function is usually defined by assuming that particles are in 
completely separated phase-space {yi 7^ 1/2) implicitly. However, the function can not be 
applicable when particles are found in the same phase-space gap {yi = 1/2) ■ We introduced 
the constant parameter b to take this into account [IB]. The correlation length can be 
obtained by the integration of correlation function C2 by assuming the constant single 
particle density, which is defined as 

„ 1^^ C2dyidy2 , 

F2-l = K2 = — — , (4) 

/ Pi{yi)pi{y2)dyidy2 

where F2 is the second order normalized factorial moment and K2 is named as the factorial 
cumulant Mathematically, K2 coincides with the inverse of NBD k parameter. In 



4 



T. Nakamura for the PHENIX Collaboration 



the case of the correlation function parametrized in Eq. Q, a relation between the NBD 
k parameter and the integrated two particle correlation function is obtained, 

1 „ , 2a^^[67]/^-l+e-^'^/^] b 

^ 2' 



(5) 



k{5r)) d-rf 
Fig. 01 shows measured NBD k values, which are corrected by the MC simulation, as 
a function of pseudo rapidity gap sizes in the Au+Au collisions at ^/snn = 200 GeV. 
The correlation lengths are extracted by the fits to Eq. by fixing the parameter a 
at 1. Fig. m indicates the extracted parameters ^ and 6 as a function of the number of 
participant nucleons. Correlation length indicates a power law behavior in the relation 
with the number of participant nucleons. 



10 = 



PHENIX Preliminary 



10' 



PHENIX Preliminary 



10-= I i i i i ' — ^ i i 1 10-' I ' ' ' i ' — ^ ' i — I 

10' 10' 

Number of participants Number of participants 

Figure 4. Extracted parameters ^ (left) and b (right) in Eq. (jSJ as a function of the 
number of participant nucleons. Solid line in the left figure is a fit with one exponent. 
Error bars indicates the fitting errors. 

4. Summary 

We measured event-by-event fluctuations for the charged particle multiplicity in Au+Au 
and Cu+Cu collisions at ^/snn = 200 GeV and 62.4 GeV. The measured charged particle 
distributions by PHENIX agree with the NBD as well as in the past experiments. NBD 
k parameters are not scaled by the average multiplicity but scaled by the number of 
participant nucleons or system size in the Au+Au collisions for the different collision 
energies. The scale dependence of NBD k parameters for the various pseudo rapidity 
gap sized are presented. The extracted two particle correlation lengths from the scale 
dependence have a power law behavior as a function of the number of participant nucleons. 

REFERENCES 

1. P. Carruthers et al, Phys. Rev. C 44 (1991) 2663. 

2. R. Korus et al, Phys. Rev. C 64 (2001) 54908. 

3. A. Mekjian, nucl-th/0503009 

4. K. Adcox et al, Nucl. lustrum. Methods Phys. Res. A 499 (2003) 469. 

5. E.A. De Wolf, I.M. Dremin, W. Kittel, Phys. Rept. 270 (1995) 1. 

6. T. Abbott et al. (E-802 Collaboration), Phys. Rev. C 52 (1995) 2663. 

7. M. E. Fisher, Physica 28 (1962) 172. 

8. See e.g. H. E. Stanley, Clarendon Press, Oxford (1971). 



